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Analytical and numerical methods are applied to principal chiral models on a two- 
dimensional lattice and their predictions are tested and compared. New techniques 



C^ , for the strong coupling expansion of SU{N) models are developed and applied to 

the evaluation of the two-point correlation function. The momentum-space lattice 
— 1 ' propagator is constructed with precision O {/3^^) and an evaluation of the correlation 

length is obtained for several different definitions. 
^ . Three-loop weak coupling contributions to the internal energy and to the lattice 

Q^ i P and 7 functions are evaluated for all A^, and the effect of adopting the "energy" 

^^ I definition of temperature is computed with the same precision. Renormalization- 

group improved predictions for the two-point Green's function in the weak coupling 
^^ ■ ( continuum ) regime are obtained and successfully compared with Monte Carlo data. 

Q^ . We find that strong coupling is predictive up to a point where asymptotic scaling 

in the energy scheme is observed. Continuum physics is insensitive to the effects of 
the large A^ phase transition occurring in the lattice model. Universality in N is 



p ' already well established for A^ > 10 and the large A^ physics is well described by a 



"hadronization" picture. 



X ■ PACS numbers: 11.15 Ha, 11.15 Pg, 75.10 Hk 



Typeset Using REVTEX 



I. INTRODUCTION 

The study of two dimensional SU{N) x SU{N) principal chiral models is strongly mo- 
tivated by the deep analogies between this class of field theories and four dimensional non 
Abelian gauge theories. We only mention here asymptotic freedom and the existence of a 
large A^ limit that can be represented as a sum over planar diagrams. A standard lattice 
version of the continuum action 

S = Jd'x^TTd^Uix)d,U\x) (1) 

is obtained by introducing a nearest-neighbour interaction 

Sl = -2NpY.^eTT [UnUU^] , /? = -i- . (2) 

n,ij, 

No exact solution of these models is known, even in the large A^ limit. An exact S-matrix 
has however been conjectured 0,@,0, and numerical evidence seems to indicate that the 
corresponding bound state spectrum is reproduced in the continuum limit Q]. Also the 
mass-A parameter ratio has been conjectured by using the Bethe Ansatz approach, and the 
result is [|1 

_ M _ /Stt sin tt/N 
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Information coming from S-matrix and large A^ factorization leads to the conclusion that 
when N ^ oo principal chiral models are just free field theory in disguise. In other words, 
a local nonlinear mapping should exist between the Lagrangian fields U and some Gaussian 
variables 0. The nontriviality of the realization may however be appreciated when con- 
sidering the two-point Green's function of the Lagrangian field: while at small Euclidean 
momenta there is substantial evidence for an essentially Gaussian (free field) behavior, at 
large momenta, where results from standard weak coupling perturbation theory are expected 
to hold by asymptotic freedom, the Lagrangian field seems to behave more like a composite 
object formed by two elementary (Gaussian) excitations, which cannot however appear like 
deconfined free particles. This elementary "hadronization" picture, where however the La- 
grangian fields themselves play the role of hadrons, is strongly supported by all numerical 
evidence that we have produced, and for A^ ^ 6 is universal, i.e. independent of A^, which 
confirms that the 1/A^ expansion, were it available, would be an extremely predictive tool 
in the analysis of these models. 

In the persistent absence of such an expansion, we here tried to apply all other analytical 
and numerical methods of lattice field theory that seemed appropriate to the problem at 
hand. In particular we sistematically extended the strong coupling character expansion to 
SU{N) groups, generating compact formulas for many "standard" integrals and character 
coefficients and extending the series for the free and internal energy, the magnetic suscepti- 
bility and various definitions of the mass gap. Since principal chiral models possess a wide 
scaling window in the strong coupling domain, these results become a powerful tool for the 



analytic computation of physical (continuum) quantities with very small systematic error 
due to scaling violations. 

We also performed many three-loop weak coupling lattice computations in order to im- 
prove our understanding of the approach to asymptotic scaling. We found that introducing 
a new coupling Te proportional to the energy [|^, the so called "energy" variable, allows an 
amazingly impressive improvement in the convergence of lattice results towards the contin- 
uum "asymptotic scaling" predictions, not only when considering such ratios as M/Al, but 
also when parametrizing the running coupling dependence of the two-point Green's function 
in the large momentum regime. 

Some of the numerical and analytical results presented here were announced (and pre- 
sented without proof) in our Refs. [@,|1. The study of the two-point Green's function in 
turn is fully original and constitutes another impressive test of scaling and universality in 
the context of large N principal chiral models. 

This paper is organized as follows. 

In Section |I| we describe a new technique for the strong coupling expansion of SU{N) 
models. The strong coupling series of the free-energy density is calculated up to O (/3^^) for 
N>7. 

In Section |T| we apply our new technique to the evaluation of the two-point correlation 
function. Strong coupling series of several different definitions of correlation length are 
presented. 



In Section |V| we compute the three-loop weak coupling contributions to the internal 
energy, and to the lattice /3-function and anomalous dimension of the fundamental field. The 
corresponding quantities in the energy scheme are also calculated. Continuum predictions for 
the two-point correlation function are obtained by solving the corresponding renormalization 
group equation. 

In Section ^ we present the results of Monte Carlo simulations for several large values 
of N, and compare them to our analytical (strong and weak coupling) calculations. 

II. STRONG COUPLING EXPANSION OF SU{N) MODELS 

As stated in the introduction, a renewed interest in the strong coupling expansion of 
chiral models was stimulated by the observation of precocius scaling well within the expected 
convergence radius of the strong coupling series 0] , and by the relevance of the complex (3- 
singularities of the partition function close to the real axis [^ . 

Strong coupling in matrix-valued lattice models was pioneered many years ago by sev- 
eral authors. Most studies were however addressed to the (relatively simpler) problem of 



computing U{N) group integrals [[10| , pT|p!2| , pr3[ , while not many results are available con 



cerning SU{N) integration |]14|,|T5],|T6| . Moreover, while very general compact formulas can 
sometimes be written, often these formulas require a lot of supplementary work in order to 
extract the directly relevant information. Alternatively, tables of numerical coefficients can 



be routinely generated by computer programs IjT^ , but generality of the results is completely 
lost. In the search for sufficiently general, but at the same time manageable results we tried 
to follow a pathway originally opened in Refs. [^,^,^, whose notation we shall try to 
follow as far as possible. 



We shall focus on the free-energy density of the one dimensional SU{N) x SU{N) chiral 
model, which can be reinterpreted as a generating functional for SU{N) group integral in 
the "standard form" 



Sp^g = fdU (Trf/f (TiU^y 



(4) 



where dU is the normalized Haar measure for SU{N). According to the definitions, the 
free-energy density F can be obtained from evaluating 



exp F(/3) = f dU exp [N(3 (Tt U + TtU^ 



(5) 



In turn knowledge of F allows the knowledge of the coefficients 2;(r)(/5) of the character 
expansion of the integrand: 



exp 



'Nf3 (Tr U + TtU^ 



exp [-F{(3)]J2d{r)Z(r)X{r){U) 
(r) 



(6) 



where J2(r) is a sum over all finite dimensional irreducible representations of the group, X{r) 
and d(^r) are the corresponding characters and dimensions. Let's now recall from Ref. |T^ a 
few exact results concerning the U{N) integration. If we denote by 



^m,N = det [Im+i~j[Sj\ , 

where Ii are the modified Bessel functions and s = 2Nf3, we find that 

Fn{(3) = lnAo,7v(s) , 



A^mAP) ^ {[detUD 



A. 



m,N 



and the following non linear ordinary differential equations are satisfied: 

2 T>r2' 



I d d 1 

s ds ds ' 1 — A^ 



N 



/rfAi, 



N 



N' 



\ ds 



Ai,iv + 1 



^1,n) ^l,N 



(7) 



(8) 



(9) 



(10) 



— (lnF^-lnF^_i) 
ds 



Ai,;v (d^i^N N ' 



l-A?^ V ds 



Moreover [^ 



a: 



m,N 



^m+l,N^m-l,N — ^m,N-l^m,N+l [^ — '^l 



,N 



(11) 



(12) 



results [1T|J20 



As a consequence of these equations, one may obtain the following strong coupling U{N) 

Ai,^(/3) = M2N(3) + O (/?=^^+2) , (13) 



A2,iv(/5) = JN{2Npf - Jm-i{2NP)Jn+i{2NP) + O (/3^^ 



(14) 



(JXJ 



4Af+4 



(15) 



fc=i 



The passage to SU{N) is effected by the relationship 

oo 

FNiP) = MP) +ln E ((dett/)') 

fc=— oo 



(16) 



implying 

oo 

F^iP) = N^P^ + 2M2NP) - 2Jr,-i{2NP)JM+i{2NP) - Y.kJl^^{2NP) + O (/?=^^) . 

(17) 



k=l 



We can now make use of known relationships concerning the series expansion of a product 
of Bessel functions to obtain a wide number of "standard" integrals: 



5, 



p,N+p 



Sp,p = P^- , P< N ; 






p< N+1 



(18) 



(19) 



Sn+p,n+p = {N + py. + Y, 



^1 (p-qV- 



{N + p)\ 
{N + q)\ 



2N + 2q-2 
q-1 



p<N+l ; (20) 



{2N + 2q)\{2N + p)\ 



''''''-' h (2^ + qnN + q + l)!(iV + q)\\q , 



p< N +1 



(21) 



One may check that essentially all Hansen's results ||T3 are correctly reproduced and ex- 
tended to arbitrary A^ by the above formulas. Eqs. ( p!8| - pl| ) are our first set of compact 
results concerning SU{N) integration. 

Even more important is however the possibility of extracting similar results for the coef- 
ficients of the character expansion. Let's briefly recall from Ref. |T^ the following expression 
of U{N) character coefficients Z(^r){P) 



C?Ai...Aiv2:Ai...Aiv 



det[JA,+,-.(2iV/j)] _ 
det [Ij-^{2NP)] 

Its SU{N) counterpart can be represented in the form 



(xUiu)) 



Z^r) 



Eoo ~ 



n 



, ^(s-lJV) 



(22) 



(23) 



where by definition 



^(s-i^) = ( (detuy ) 



(24) 



d'{r)Z(r+s-l'^) 



{{detUrxhiU)) 



(25) 



and all quantities are computed with the U{N) measure. These results become however 
rapidly useless with growing N, due to the intractability of the determinants. We may 
however easily recover from the previous results the coefficient of the fundamental character: 



NzM 



d 



d{2N(3) 



Fn{P) 



CXJ 

NP + Jn-i{2N[3) - Jn+i{2N[3) - ^ J^+fc(2iV/3) J^+fc+i(2iV/?) + O [(3^''-^) . (26) 



fc=0 



Moreover we may want restrict our attention to the first ~ 2N orders of the strong coupling 
series, in order to exploit the techniques discussed in the first part of this section. In this 
case the character coefficients of U{N) can be reconstructed by the use of Schwinger- Dyson 
equations [|1^ . Appendix A is devoted to a presentation of our main results concerning the 



compact evaluation of U{N) and SU{N) character coefficients. 

For what concerns the physical problem of evaluating the d-dimensional free-energy den- 
sity in SU{N) X SU{N) chiral models within the strong coupling expansion, we address the 
reader to Ref. |T^ for a discussion of the relationship between the SU{N) and the U{N) 
character expansions. Suffice it to say that knowing the U{N) expansion to O (/J^™") allows 
an immediate identification with the corresponding SU{N) expansion when N > m. When 
N = m only a minor modification is required in order to avoid a double-counting due to the 
self-duality of the representation (l^/^;0) = (0; 1^/^). When N < m the procedure must 
be much more careful, but in any case this condition would violate the already mentioned 
restriction on the applicability of our approach, which does not trivially extend to more than 
2A^ orders of the strong coupling series. 

In order to exhibit some new results concerning the strong coupling series for the free- 
energy of the principal chiral models, let's first introduce a notion of "potential" intimately 
related to that presented by Green and Samuel in Ref. 
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St 



For comparison, we mention that 



2^-tf^2t^2s 



N' 



2-2t 



(27) 



(28) 



where V and V have been defined in Ref. [|T8|. The tenth order strong coupling character 



expansion for the free-energy in d-dimensions is then pT 

d 



T= r^ + 2 



[z{ + 1^42) + 



A\VA\ 



'A^ 



96 [3] 1^21^? + 



24 



Vnzl + W^izt 


+ 




^ + 576 ^ + 1536 ^ 


\ / \ "^ 


/ 


\ ^/ J 



^Ul296 f 



^1^ + 



Wu zf + 



56(2 ) +4656(3 )+46272(^ 



95232 



6 



ao 



(29) 



Eq. (H) holds for all U{N) and for SU{N) when A^ > 5. The SU{N) result for the 
"potential" Wst can be obtained in the form 



2* 



1 + ^JnC^N/S) 



p 



2* 



Jn+2{2N(3) ( 1 + 



N^ 



-^ \ t-s 



1 

N 



-^ \ t-s 

l.t- 

N 



211-]^^*"-^ 



+ 



+ JN-2i2NP) (1 - -) - 2M2NP) (1 - — ) 



t-s 



(30) 



By explicitly expanding in powers of P the 14*'^ order character expansion for the free energy 
(presented in Ref. ^ up to O {p^'^) ), we obtain for U{N) models (A^ > 7) 



+ 
+ 



2/52 + 2(3^ + A/3^ + 



14 + 



N^{5N^ - 2) 

(iV2 _ 1)2 



/3' + 



56 + 



]N^{5N^ - 2) 

(Ar2 _ 1)2 



/3 



10 



8A^2(^35A^2 _ 17) 2Ar2(14Ar6-llA^4^8A^2_2) IQN^{9N^ - 26N^ + 8) 

248 H 7-—^ —-: \ ——: 7— I h 



1176 + 



(Ar2-i)2 (iV2-l)4 3(iV2 - l)2(iV2 - 4)2 

432A^2 ^2N^ 2A0N^{5N^ - 2) iV2(248iV4 _ iun'^ + 43) 



/3 



12 



iV2 - 1 (Ar2 - 1)2 



{N^ 



A^2(436Ar6 - 344A^^ + 208iV2 - 48) GAN^dN^ - 26N^ + 8 



(Ar2 - 1)4 
and for the SU{N) models (A^ > 7) 



T = T ^ ^4^/3^ + 



(A^2 _ 1)2(]V2 _ 4) 



(Ar2 - 1)3 



16 



(31) 



2A^ 



Ar+2 



AT! 

^(A^ + 2)A^ 



4Ar^ SAT^-i 

+ 



TV 



+ 



(A^ + 1)! A^! 

4JY7V+1 24A^^-i 



/3^+2 + 



(A^ + 2)! (A^+1)! (A^-l)A^! 



+ 



A^! 



+ 



4A^ 



N-2 



{N -2)\ 



/? 



Ar+4 



+ 



2Ar^+4 4Ar^+2 (iv + 4) 4ivJV+i ( jv2 + 4]v ^ 2) 8Ar^+i(l + A^ - Ar2 



3(A^ + 3)! (A^ + 2)! (A^ + 1)(A^+1) 

4JYJV+3 24A^^(A^ + 2) 4A^^ 



(Ar2_i)(jv + l)! 

\N^ 24Ar^+i 
+ 



+ 



3(A^-1)(A^-2)A^! (A^ + 1)! (A^-l)! A^! (A^ - l)Ar! 

4A^^ 112A^^-i 24A^^-2- 



^A^+6 ^ Q |^^A.+8^ _ 



(A^-l)2(A^-3)! A^! (A^-2)! 

The S't/(6) result, including an analysis of the O ii?^) contribution, is finally 

o ^2 o ^4 56 . 84038 ^„ 459308 ^^ 548436429 ^.o 
^ M ^ g /^ ^ ^225 '^ 1225 '^ 85750 '^ 



(32) 



(33) 



The strong coupling series of the internal energy E can be obtained by 



E 



1 - 



Adf3 



(34) 



In order to extend the strong coupling series of thermodynamical functions to higher 
orders and in order to compute the strong coupling series of correlation functions, it will prove 
convenient to define more involved potentials corresponding to nontrivial loop topologies 
and three-body interactions. For the purposes of the present paper, and in view of further 
developments, we shall define 



W, 



ab 



a+b I A^2 
(1;1) ^ 



'(i;i)' 



d{i-i) 



-{2;0) 



d, 



+ 



(2;0) 



^(12;0) 



+ z: 



dn 



(i;i)" 



(i;i) 



^(2;0) 

Am 



+ 



^(12;0), 



2zf 



2a+2b 



(35) 



Wabc = ^(l;l) (^{2;0)'^(2;0) + ^(l2;0)'^(l2;0)) (^('2;0)'^(2;0) + ^ll^;0)d{l2;0)) - 2N'^Walzl 



2b+2c 



- Z, 



^b+c 



(i;i) (-^(2^0) '^(2;0) + '2(i2':o)'^{i2;0) j + permutatious of (a, 6, c) 



+ ^(1;1) V"(l,l) ~ "{i;i)J - 4iV Z-^ 

and the three-body potentials 



4 ^2a+2b+2c 



K 



qst 



N 



l-2t 



{^h-fiAm + ^^2,l;0)'^(2,l;0) + 42,1)'^(2;1)) ^f2;0)'^(2;0) 

+ (^^2,1;0)'^(2,1;0) + ^('i3;0)'^(l3;0) + ^^i2;i)'^(l2;l) j ^fl2;0)'^(l2;0) 
(^(2;1)'^(2;1) + ^('l2;l)'^(l2;l)j ^fl;l)'^(l;l) 



+2 



1-t 



2iV2*+i4W^< 



g+s,i+l 



,2-t^2t+1^5+2s^^^^ 



22-t ^T2t+3 3q+2s 



(36) 



(37) 



III. THE STRONG COUPLING EXPANSION OF CORRELATION FUNCTIONS 

A typical application of the strong coupling analysis amounts to the evaluation of the 
so-called "true mass gap" of the model, which in turn is defined to be the coefficient of the 
asymptotic exponential decay of the two point correlation function of the order parameter. 
In SU{N) X SU{N) chiral models one defines 



G{x) = l{TT[U\x)Um) 



and the true mass gap is 



fj, = — lim 



In G{x) 



|x|^oo \X\ 



(38) 



(39) 



It is also possible to introduce the lattice moinentuiii transform 



G{p) = Y.G{x)expU^^^) (40) 



and extract the mass gap from the (imaginary momentum) pole singularity of G{p), i.e. by 
solving the equation 

G-\p = ifi) = . (41) 

There are a few well-known (but not necessarily well understood) results concerning the 
evaluation of the mass gap. In particular it is often stated that in the strong coupling 
regime G{x) "does not exponentiate" and it is therefore necessary to define the wall-wall 
correlation 

G^(x|i) = ^G(x^,X||) (42) 

x± 

enjoying the exponentiation property, and assume 

/. = - hm i^^^4^ . (43) 

|x|||^oo |a;||| 

Appendix B is devoted to a short discussion of this question, which finds an easy illustration 
in the context of the exactly solvable Gaussian model. 

In practice, when extracting physical quantities from Monte Carlo simulations, one typ- 
ically faces a situation where it is non- realistic to solve numerically Eq. (^T]), which requires 
an analytic continuation to negative p^ of the (usually poorly known) function G~^{p). We 
shall be able to perform this exercise in our specific study, but this is far from being an eas- 
ily generalizable practice. On the other side, working on a finite lattice often prevents from 
exploring a sizeable region where exponentiation may hold with small errors for arbitrary 
/?. We therefore found convenient to analyze explicitly the strong coupling expansion of the 
Green function G{x) for finite fixed lattice distance, in order to establish a benchmark both 
for exponentiation and for the analysis of the numerical results. 

Without entering many details of the (sometimes cumbersome) strong coupling calcula- 
tions, we want to sketch the main ingredients and logical steps before presenting our results. 
In passing we notice that we might consider a slightly more general Green's function than 
Eq. (HD: 

{Xir)[UKx)Um) u,. 

G(r) = ^ , (44) 

"(r) 

where r can be restricted, for principal chiral models, to the completely antisymmetric 
representations (1";0). The strong coupling expansion is naturally ordered in the length 
of the path that we choose in order to connect the point x = (xi,X2) with the origin. In 
turn, since every path can be decomposed into a sum of shorter paths, we may in general 
establish recursive relationships connecting the (fixed length) contributions to G{xi,X2) with 
the (shorter length) contributions to Green's functions of nearby points. 

9 



The strong coupling character expansion of Green's functions involves, to lowest orders, 
only a summation over properly weighted self-avoiding walks. The effect of "potentials" is a 
higher order contribution that must be included by considering bifurcating paths "dressed" 
with properly chosen representations of the link operators. For all U{N) groups (A^ > 2) and 
SU{N) groups (A^ > 4) the O {(3^) strong coupling character expansion may be represented 
by 



G{XI,X2) =^^+"2 



Co{xi,X2) + C2{Xi,X2)zl + Ci{Xi,X2)zl 



+ A{xi,X2)zlWii + B{xi,X2)W2i + 0[zl 



(45) 



where the quantities C2k{xi-,X2) represent the number of self-avoiding walks connecting the 
origin with the lattice point x = (a;i,a;2) and whose length is / = Xi + a;2 + 2k, while the 
functions A and B satisfy the relationship 



A{xi,X2) + 2B{xi,X2) = 2{xi+ X2)Co{xi,X2) 

In Appendix C we derived the following results 

Co{Xi,X2) = 





(46) 



(47) 



C2ixi,X2) 




a^i(a^i + 1) 3^2(3:2 + 1' 
a;2 -F 1 xi + 1 



(48) 



C4 (3:1,3:2) 




X1X2 + 2xi + 2x2 + 
l)(a;2 + 



2(xi + l)(xi + 2) 



+ 



X2(x2 + 3) Xi(Xi + 3) 

Xi + 1 X2 + 1 

(xi -l)xi(xi + l)(xi + 2) 
2(x2 + l)(x2 + 2) 



2xiX2 
Xi + X2 



, (49) 



A{xi,X2) 



,Xx "T X2 
Xi + X2 



Co(xi,X2) 



(50) 



E(xi,X2) 



2XiX2 
Xi +X2 



Co(Xi,X2) 



(51) 



Writing down and solving recursion equations for the coefficients of the strong-coupling 
series for Green's functions is a very powerful but by no means simple or efficient approach 
to the evaluation of the mass gap. In practice we may observe that, to any given order of the 
strong coupling expansion, the recursive relations between coefficients imply that only a finite 
number of short distance wall-wall correlations for L < L may violate exact exponentiation 
because of boundary condition effects. Therefore one may compute the quantities 



10 



lnG'^(L + l) - \nG,,{L) (52) 

for the first few values of L until they become constant. The constants obtained are directly 
related to the masses by the relationships 

Hside = \nGside{L) - InGsideiL + l) L>Ls , (53) 

^ = \nGd^ag{L) - \nGdiag{L + 1) L > L, , (54) 

where by definition 

oo 

GsideiL) = E GiL,X2) , (55) 



Gdiag{L) = Yl G{L -X2,X2) ■ (56) 

X2=—oo 

It is possible to prove that the alternative definition of /x based on solving the equation 

G-^{p = i^) = (57) 

perturbatively in the strong coupling expansion parameter gives stable results only after all 
correlations up to L = L have been included in the Fourier transform of the propagator, and 
the result obviously coincides with Eqs. (0) and (|5^) . It is therefore convenient to construct 
explicitly the Fourier transform of the inverse propagator and extract the relevant physical 
parameters by a direct analysis of this last quantity, making use of the above considerations 
in order to establish the accuracy of the computations, which in general will not coincide 
with the precision reached in the evaluation of G~^{p). 

The evaluation of G~^{p) is dramatically simplified by the observation that any strong- 
coupling expanded two-point Green's function can be unambiguously separated into the 
form: 

G{x) = Go{x) + AG{x) , (58) 

where Go{x) is originated by the Fourier transform of the generalized Gaussian propagator 

'^''^^ ^ Ao{z,) - 2zi5o(zi)E^cosp^ ' ^^^^ 

and the functions ^0(2:1) and -Bo(^i) are uniquely determined by the conditions 

Go{0) = 1 , 
Go{m = l-E{z,) = \^ ^ e{z,) , (60) 

which can be cast into the form: 

11 



1 = ^vi^2ny /25o\'" 

zB, = ^ . (61) 

From the results of the previous section we may extract: 

e{zi) =zi + 2z\ + ^z\ + 2z\Wxi + ^1%z\ + 12^1^11 + 2^11^31 

+ 140;z? + l^zlWxx + e^i^VTso + S^^VTsi + 2^131 + ^iV^ioo + O [z'A . (62) 



As a consequence we obtain 

Ao(zi) = 1 + 4^1 + 124 + 60^1^ + \<oz\Wrx + +316^1^ + 96^1^1^11 + \<ozlW'>,r + 1844^^° 

+848^^W^n + ^^z\Wl^ + 48^1^1^20 + 64ztW^3i + X^z^Vxzx + 84\/ioo + O izf) , (63) 



+6z^iyi'i + 6z^W^2o + 2zi'iy3i + 2^1-^131 + ^iVioo + O izf) . (64) 



A direct evaluation of all coordinate-space Green's functions that are non-trivial to O [z^^ 
allows us to determine AG(x) with the same precision, since Gq{x) is easily obtained by 
antitransforming Eq. ([59|) . By noticing that the first non-trivial contributions to I\G[x) are 
O (2^), it is now relatively easy to evaluate directly G~^(p) to O {z^\ We obtained: 

G-\v) =A{z^) + z^B{z^)f + z\C{z^)v\v\ 

+ z\D{z,) (v\ + V% + z\E{z,)fv\vl + O (4^) , (65) 

where v\ = 4sin2 (P/./2), p^ = E^P^, and 

A(zi) = Ao(zi) - 42i5o(2i) - 8A(^)2i^ (1 - 4^1 + 3^1^ + ^z\ - 3^^) 

-16A(2)zi« (1 - 4zi + 2z?) - 8A(3)^9 - 8A(^)^f - 16A(5)z10 + O (z"^) , (66) 



B{z^) = Bo{zi) + AA^^'^zl (1 - 2zi + 2z\ + -izl) + 8A(^)4 (1 - 2zi) 

+10A(3)^8 + 8AWz;' + 8A(5)2;' + O (^1°) , (67) 



C{zr) = -2A(i) - 4A(2)z2 _ sA^^)zl - AA'^'^h^ + O [zl) , (68) 

D{zi) = -2AW - 2A(=')^i - 2A(4)^2 ^ O (^3) , (69) 

E{zi) = A(3) + O [zl] . (70) 
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We have defined the following combinations of potentials: 

A(i) = ^-4^/21 - 2zfWii - 1 , (71) 

A(2) = _z:^f^W:ii + 4:Z^^W2i + z^'^WiQ-7z^^Wii-4: (72) 

A(=^) = z:^'^Wii + 2z:^^W-ii~Az^^W2i-2z^^Ww + Az^^Wii + l (73) 

A(^) = 2z^^W2i + z^^Wiii - z^^Wii - Az^'^W^i + z^^Wai + lQz^^W2i 

-2z^^W2o - "iz^^Wl^ - 2QzY^Wu + ^z^^Wio - 9 , (74) 

A(5) = z^'Vuo + ^^r'V22i - zfv.si - ^r Vioo - ^r'H^ii + z^^w^^ - s^r'w^si 

^5 q c^Q 

+ y^r'W^2i - 3^r'W^20 - 2^r'W^n - SS^r^^n + IS^r'W^io - y • (75) 

Eq. (|65|) is a rather compact collection of physically relevant results concerning the strong 
coupling regime of the models under investigations. We notice that the nearest-neighbor 
Gaussian structure of the propagator starts being violated to O (zf ) in diagonal correlations 
and only to O (2^) in side correlations. We therefore expect a substantial agreement of the 
ratio ic/iw, where ^q = {x^) and ^w = ^/fJ^side, with its Gaussian value even for not too 
strong values of the coupling (as we will see in Section |Vp. 

By standard arguments we can establish relationships between the propagator, the sus- 
ceptibility and its second moment: 

00 

X= E G'(Xi,X2) = Gs^deiO) + 2 ^ G',,,e(^) 
3::i,a;'2 L=l 

oo -| 

= G,..,(0) + 2 X: G,.aAL) = j^ , (76) 



1 oo 

X {x'^) = J J2(^l + xl)G{xi,X2) = J2 L'^GsideiL) 

a;i,X2 L=l 

Hence we obtain 

M'g ^ t\t = ^^ , (78) 

^ (x2) z,B{z^) ^ ' 



'"^ ZiB{z,) 



13 



(79) 



Moreover by solving appropriate algebraic equations we find 

ML. = 2 (ch/.... - 1) ^ M^ + "'"^^"^l^}"'^ + ... (80) 

7W,,,^ - 4 l^ch-^ - IJ _ 7W^ + -^^ l^^2~ + ~^) 45(^0^ + - 

(81) 
By substituting our explicit results we find 

X = 1 + 4^1 + 12^2 ^ 3g^3 ^ ^QQ^4 ^ ^284 + 82i-=^Vrii)^i^ + (788 + 4821-^1^11 + 8A(^) 



+ (2204 + 216^1-2^^11 + 8^1-^1^31 + 32A(^))4 + (6068 + 800^1-21^11 + 48zi-^1^3i + 88A(^) 
+ 16A(2))28 + (16820 + 8zi-Vi3i + 4^1-^1/100 + 290421-^1^11 + 24^1-^1^21 + 24^1-^1^20 
+200^1-^11^31 + 304A(^) + 64A(2) + 8A(^))2?+(46172 + 482rVi3i + 242i-Vioo + 984021-^11^11 
+ 192^1-^11^21 + 144^1-^11^20 + 70421-^11^31 + 1000 A^^) + 160 A(2) + 64A(=') + 8A(^) + 16 A^^))^!^ 

+0 (zl') , (82) 



2iM| = 1 - 421 + 322 + (2 - 221-211^11) 2^ - 4A(i)2i^ + (6 + 621-211^11 - 221-^11^31 + 16A(i))2? 

-(l6A(i) + 8A(2))4 + (14 - 22i-Vi3i - 21-^1/100 - 42i-2ll^ii - 221-^11^21 - 621-^11^20 
+ 10^1-^11^31 + 4A(^) + 32A(2) - 10A(3))2i^ + (12A(^) - 16A(2) + 32A(3) - 8A(^) 
-8A(^) + 162r2ll^iiA(i))2i' + O (zl^^) , (83) 



2iZg = 1 - 22 + (-6 - 221-21^11)2^ - 4A(^)2i^ + (-18 - 22i-2ll^ii - 221-^11^31 + 8A(^))2^ 
-8A(2)4 + (-98 - 22i-Vi3i - 21-Vioo - 5221-211^11 - 221-^11^21 - 621-^11^20 + 22r^ll^3i 
-28A(^) + 16A(2) - 10A(3))^8 + (60A(i) + 16A(2) - SA^^) - 8A(5) + 1621-21^11 A(i))2;' 
+0 (zl') , (84) 



22 

^is^de = - In 21 - 221 " (2 + 221-21^11)2^ ' (^ ^ ^^i^^^l + 4A(^))2i^ 

-(4 + 221-^11^31 - 6A(i))26 - (y + 421-211^11 + 421-^11^31 - 16A(i) + 8A(2) + 2A(3))2[ 
-(28 + 3221-211^11 + 421-^1^21 + 621-^11^20 - 421-^11^31 + 22rVi3i + 2rVioo + +28 A^^^ 
-16A(2) - 2A(3) + 2AW)2^ - 2A(6)2^ + O (2?) , (85) 
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/^'i-3 _ _ i^2^i _ zl - (I + 2zfW,, + iA(i) 



V2 ^ ^ V2 



/25 



^ g + 4^r2^n + 2zl'W,, - 2A(i) + A^^) + Ia^^))^^ + Q {zl) . (86) 



The contribution proportional to 

A(6) = z-^Wsi - 2z\W2i + ^ilVii (87) 

in Eq. (|85|) was obtained by evaluating the O (^jf^) contribution to (7(3, 0) and could not 
have been predicted from the knowledge of the O i^z^) contributions to G~^(p): this is in 
accord with our previously discussed considerations. Further useful results are: 

\^Gs^de{^ = Llnzi + 2{L + l)zi + -{L + 7)zf + 2{L + 2)^ + -{L + 51)zl + 2LzlWn 

6 5 

-4(L - l)zlWu + 4:LziW2i + 2(5L + 9)^^ + 4(3L + l)z^Wu - 2(3L + 1)^12^21 + 2LW31 
+2zfA^^^SL,o + ^(172 - 3L)zl + 4(L + 1)2^1^10 - 12(L - 1)^1^1^11 + 4(2L + 3)zlW2i 



+2(L - l)zlWii + 2zI{A^'^^ - 2A(^))5l,o + 22(2L + 5)zl - 4(L + 4)^1^1^10 + 2(15L + 28)^1^1^11 

-2(L - 3)4W^n + 2(L + 2)ztW2o - 20z^W^2i + 2{L + Q)zfW3i + 2{L - 1)1^41 - 2(2L - 3)z^W^ii 
+4(L - 1)^211^21 + 2(L - l)zfWui + 2z^^LV3i + zfLV^oo + 2zfA^^^SL,i 
+2;2^(2A(6) + AW - 2A(3))5l,o + O " 



25 1 

\nGd^agiL) = L lu 2zi + (L + A)zf + (2L + —)zt + (L + l)zl Wn + ^{L - 1) W21 

+^ztA^'^6L,o + ^-^^^4 + ^4Wio + (2L + 9)41^11 + (L + 4)^1^21 + hw^i 
+ ^z!mi + U^^'^Sl,! + U(A^'^ + 2A(2) - 4A(^))5z.,o + O (zf) . (89) 



The results we have obtained are expressed in the language of the character expansion. 
However within the precision of our expansion it is possible, for all U{N) and SU{N) groups 
with A^ sufficiently large, to convert the results into a standard strong coupling series by the 
following replacements (dictated by our previous analysis) 

Zl ^ P , (90) 

^r'w^ii ^ , (91) 

47V^ 
z^^W,o - 2z^^W2i ^ j^^—^ ^-^^4 , (92) 

2N'^(5N'^ — 2) 
zi'W2o ^ 2z-,'W,i ^ (^2 _ 1). ,— 10 ' (93) 
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-8w ., NH9N'-eN' + 2) 



-1 M'" - (A,. _ 1)3 .~J • (''4) 



_2~ _ N^{7N'^-2) 



-. "^"--(1^131)7^.-7. (95) 

.— N^(llN^-QN'^ + 2) , , 

'""'- - (TV. -1)3 ^ .-ZL 11 • (''«) 

_4Ar4 

^fVioo ^ , (98) 

%-Vno =^ , (99) 

%-Vi3i ^ , (100) 

In Appendix D we give explicitly the N = oo strong coupling series of some relevant quan- 
tities considered in this section. 



IV. WEAK COUPLING EXPANSION 

For both continuum and lattice, at low temperature, the perturbative expansion is per- 
formed by setting 

U = e'^ , A = Y,TaAa , (102) 

a 

{Ta are the generators of the SU{N) group and Aa are A^^ — 1 real fields) and expanding U in 
powers of A. The above change of variables introduces an additional ill-defined determinant 
in the partition function, indeed 



[dU] = K X{[dAa]e-^- , 

a 
Q ^Y-T^ 1 2(1- COS i„) 

Sm = -7T Z^ Tr In J- , (103) 

^ n. At 
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where Abe = Y.a'^'fbacA.a and K is an irrelevant constant. In dimensional regularization, the 
measure term does not contribute, as a consequence of the rule / d'^k = <S=^ 6'^{0) = 0, 
where d is the space dimension. 

Short weak coupling series for the free-energy density of U{N) and SU{N) chiral models 



on the lattice were presented in Ref . p5[ . We calculated the energy density up to three loops 
finding 



E 



where 



1 - ( -Re Tr [f/.t/I+J ) 



N^ -2 



N' 



IN^I3 



1 + 7j + ^ + - 



(104) 



Oi 



3Ar4 - 14A^2 ^ 20 N^-AN^ + 12^ N^-SN^ + 2A^ 



768A^4 g4]v4 

Qi and Q2 being numerical constants: 



64A^4 



(105) 



Qi 



jPk_ 

(2^ 



k^ 



d^p 
(27r)2 \ g2^^^ P 



2 

2 ~ ^2 



Q2=E 



fi,i/ 



5{Y^h 



{ki + k2)Ak^ + k:,\ 



h.2 i,2 1,2 i,2 
A/i rLO ''-'S '^A 



- 2 



n2 



A;2 



(106) 



Qi = 0.0958876 and Q2 = -0.0670. 
We calculated 



G{T,x) 



J[dU] iReTr[f/(0)f/(x)t] exp(-5/T) 
IldU] expi-S/T) 



(107) 



in perturbation theory, and in two regularization schemes: dimensional and lattice regular- 
izations. In dimensional regularization and in the x-space 



Gd(T,x, e) 



N^ -1 T 



1 



iV2 _ 1 iv2 - 2 T2 



2N {27rYSd ex' 2N 8N {27iy<^Sj e^x^' 

1 



+ 



(108) 



where Sd = [{271^/^2' /^T {I + e/2)J and e = rf - 2. In the p-space 
N"^ -IT 



GD{T,p,e) 



2N p2 



N'^-2^^p' .^ooJ-N^ + 21 iV^l 3N^\ 



(109) 



On the lattice we obtained (neglecting O (a^) terms) 



GLiT,x,a) = 1 + -^^TE{x/a) + -^^ ^^T'E{x/a)[l + 2E{x/a)) + ... 



27V 



2N 16N 



(110) 
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where 



1 / 3 \ 

E{x/a) = — Mna/x-7£;- -ln2J 



111) 



GL{T,p,a) = 
' -N^ + 2 



T 



lem 



N^ -IT 

2N f 

B{paf + — 



N'^-2 f l\ 



4Ar2 



32A^2 



-B{pa) - -^C{pa) 



5A^4 - 25Ar2 + 30 Ar2 \ 

768iV^ + -G, +... (112) 



where 



B{pa) 
C{pa) 



W, 5, 
— m pa m 2 

27r V 2 



1 



(2 



TT 



In pa In 2 



(113) 



and Gi is a numerical constant: Gi = 0.04616363 [p6ip7| . The equivalence of the SU{2) x 
SU{2) chiral model to the 0(4) a model allows a check of these expressions, indeed for 
N = 2 they must give (and indeed they do) the same results obtained for the corresponding 



correlation functions in the standard lattice 0(4) a model [ p6|] . 

On the lattice, neglecting O (a^ In" a) terms, the correlation function Gl{T,p, a) satisfies 
the renormalization group equation 



-aS- + Pl{T)4^ + 1l{T) 



da 



dT 



GL{T,p,a) = 



(114) 



where the functions Pl(T) and jLiT) tell us how the temperature and the field U should 
vary with the lattice spacing a to keep the renormalized quantities fixed: 



PLiT) 



-a-—! 
da 



boT^ - biT^ 



h,T' + ... 



(115) 



d 



7.(T) ^ a-lnZ^ 



7iT + 72,T2 + 73,T3 + ... 

Zjj is the function entering the renormalization of the two point function 

GR{t,x,fi) = Zu(T,a^y^ GL(T,x,a) , 



(116) 



(117) 



where t is the renormalized coupling connected to the temperature by the relation T = tZt, 
and // is an energy scale. In Eqs. (|115|) and (|116|) Bq, bi and 71 are universal coefficients 
independent of the regularization scheme, and appearing also in the renormalization group 
equations giving the behavior of the renormalized quantities when varying n keeping the 
bare quantities fixed. 



bo 

7i 



N 
ANtt 



bi 



N^ 



1287r2 



(118) 



The coefficients 62^, , 72^ and 73^ can be calculated using the procedure described in Ref . 
]. We determine the renormalized functions Z^^^{T, afi) and Zif^{T, afi) that satisfy the 
equations 






;ii9) 



where t and G^ (t, x, /x) are respectively the coupling and the correlation function renor- 
malized in the MS scheme. Renormalizing in the MS scheme the expressions (|108| ) and 
(|109|) we find 



N^-1 t 






1 + .^^ ,^ .. + Oif-) 



128 (27r)2 

where 7^ is the Eulero constant. Then by imposing Eq. ( |119| ) we obtain 
zf^{T,afi) = 1 + LiT + L2T^ + 0{T^) , 



where 



Li = d (In a/i + di) , 
Ci = bo , 

5 iV^-2 

di = m 2 — TT -— , 

^2 2A^2 ' 

C2 = &1 , 



(120) 



;i2i) 



(122) 



Z^'^{T,aii) = 1 + MiT + 



W^t) "' - "^'^ 



rp2 ^ Ql^rj^S^ 



where 



(123) 



Ml = ei (In a/i + /i) , 

ei = 7i , 

The ratio of the A-parameters Aj^ and A^ is given by |jl5| 



^ = exp(-di) 

Ar, 



32exp TT- 



N^-2 ' 
2m 



(124) 



(125) 



72^ is easily obtained from Eq. ( |116[ ): 
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72i = &o7i(/i-c?i) 



{N^-l){N^-2) 



By definition, when keeping the bare quantities a and T fixed we must have 



(126) 



/i|-t = -t/i|- In Zf^(T,a/i) = /%^(t) = -b,t^-b^t^-h—t'' + ..., (127) 



and 



^'Ji:i^'^^u'{T,afi) = 7msW = 7it + 72^t2 + 73^t' + ... 



(128) 



where [3jjg{t) and 'Jjjsi't) are respectively the /3 function and the anomalous dimension of 
the field U in the MS renormalization scheme. On the other hand, deriving with respect to 
a, keeping the renormalized quantities /i and t fixed, we must obtain the (3 function and the 
U field anomalous dimension defined in Eqs. ( 115 ) and 



Pl{T) 



and 



-Ta^\nZ''^(T,afi) 
da 



d 



7l(T) = a-^lnZ/f^(T,a/x) 
da 



(129) 



(130) 



By comparing the perturbative expansions of Eqs. ( |127| ) and ( |128| ) with those of Eqs. ( |129| ) 
and ( [1301 ), the coefficients 62^ and 73^ can be written in terms of the quantities introduced in 
Eq. (|121|) and the corresponding coefficients &2^5 12— and 73 , which have been already 



calculated [p8|] 



'MS ' ' MS 



3Ar3 I 
512 (27r)3 ' 



;i3i) 



and P9| 



72 — = , 

' MS ' 

_ ?,N{N^ - 1) 1 
^3ms ~ ^ (27r)3 



(132) 



Indeed 



62^ = h2—- + 60 hi {di - do 



(133) 



and 



In particular we find 



73^ + 71 bl {di - hf - h {d2 - h) 



(134) 
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b2r 



1 N^ 



(27r)3 128 



1 + 7r- 



Ar2 



2m 



TT^ I T^. + AGl 



QN^ 



(135) 



Having calculated 62^, we can evaluate the first perturbative correction to the two loop 
relationship between the lattice scale A/, and the temperature 



A, 



-bi/bi ^ 



%T) "' "exp 



(87r/5)^/' e-S'^^ 



1 



6nT 



1 + 



1 + 



hi - &0&2, 
Nhl 





^-1 + 0(^2) 





(136) 



In the large A^ limit 



h\ - bob2^ _ 



0.0605095 . 



(137) 



In order to solve Eq. ( |114[ ), let's introduce the dimensionless function HiiT^pa) 
p'^Gl{T,p, a). A formal solution of Eq. (|114| ) is given by 



HLiT,pa) = HL{e,l) exp / dz 



e 



1l[z) 



;i38) 



where = Q{T,pa) satisfies the equation 



y^Q{T,y) = /3L(e) . 



(139) 



Defining the function 



.(r).exp(/dT|(I|) 



J^-'Ti/bo 



1 + 



71^1 - 72i&o 



/)2 



T 



+ 



l2Lhbl + 7i(^2^^o - bl)bo - 73^ &g + (7161 - 72^ &o) 

2hl 



■T^ + O T 



HL{T,pa) can be rewritten in the form 

HL{T,pa) = z{T)-^ z{Q)Hl{Q,1) . 
Solving perturbatively Eq. ( |139|) we find 
1 



e 



bou 



hihiu h\\\iu{]xiu — 1) + h2,^hQ — h\ (hfu 



hlu 



bt 



u^ 



w 



(140) 



:i4ii 



(142) 



where u = ln(p/Ai). The perturbative expansion of Hl{Q,1) can be obtained from Eq. 
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In order to get a more accurate description of the approach to asymptotic scahng we 
performed the change of variables suggested by Parisi 0, defining a new temperature Te 
proportional to the energy: 



Tf 



8N 



iV2-l 



E 



I3e 



NT, 



(143) 



The ratio of A^;, the A parameter of the j3e scheme, and Kl is easily obtained from the two 
loop term of the energy density: 



Ke ( N^-2' 

AT = ^"P V^m- 



(144) 



Within the jSe scheme, the two point function G^{Te,p, a) must satisfy the renormalization 
group equation (neglecting O (a^ In" a) terms) 



'"i + ''-'^-'^ + ^'^<''^' 



GE(TE,p,a) = 



(145) 



The /3-function of the Pe scheme can be written in the form 

rl S A/^2 



where 



CiT) = If^ 

^ ' NdT 



(146) 



(147) 



is the specific heat and T must be considered as a function of Te- Expanding perturbatively 
Eq. ( |146| ) and using Eq. (|104|) one finds 



&2b = h^ + N'^bo (a2 - a?) + Nhai 
The scale A^; is related to the variable Pe by 

, bj-boh 



Ae = (87r(3Ef^ e~^^^^ 



Nm 



''^Pe' + O 



;i48) 



(149) 



In the large A^ limit 



bj - bob2^ ^ 



-0.04009 



The function 'Je{Te) is easily obtained from the relation 

7e{Te) = 1l{T) , 
and therefore 



(150) 



(151) 
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{N^-l){N^-2) 

72. = 72. - 7iiVa, = ^^^-^^ , 

73^ = 73^ - 2-f2,Nai + -fiN^ (2a\ - as) . (152) 



The renormalization equation ( |145| ) can be solved following the same procedure used in 
the standard scheme. Defining HEiTE,pa) = p'^GE(TE,p,a) we obtain 

HE{TE,pa) = ze{TeY^ ze{Qe)H{QeA) , (153) 

where 



ze{T) = exp I / dTi 




(154) 



and Qe = ^EiTE.pa) satisfies the equation 



y^QE{TE,y) = I3e{Qe) , (155) 

whose perturbative solution is obtained from Eq. ( p.42|) by substituting 62. ~^ ^2^ and A^ -^ 
Ae- The perturbative expansion of He{Qe, 1) can be found by reexpressing perturbatively 
T in terms of Te in the r.h.s. of Eq. ( p,12| ). 

Notice that having performed only a coupling redefinition T -^ Te, it must be 
Ge{Te,P,(i) = GL{T,p,a). But when evaluating them at some finite order in perturba- 
tion theory we get different results, which are different approximations of the same quantity. 



V. NUMERICAL RESULTS 

In order to investigate numerically the large A^ limit of SU{N) x SU{N) chiral models, 
we performed Monte Carlo simulations for several large values of A^. We will show numerical 
results for A^ = 6, 9, 15, 21, 30. In Ref. [^ some results at A^ = 6, 9, 15 were already presented, 
and in the following we will use some of those data. A summary of our new large statistics 
Monte Carlo results is presented in Table |. 

In our simulations we used local algorithms containing overrelaxation procedures. We 
implemented the Cabibbo-Marinari algorithm |3^ to upgrade SU{N) matrices by updating 
their SU{2) subgroups. In most cases, the SU{2) updatings were performed by using the 
over-heat-bath algorithm |3l[] (for the "heat bath" part of it we used the Kennedy-Pendleton 
algorithm |3^). A sweep consisted in updating a number of SU{2) subgroups at all sites 
of the lattice. For relatively small values of A^ (A^ ^ 6 say) we chose to update the A^ — 1 
diagonal subsequent SU{2) subgroups of each SU{N) matrix variable. At larger A^ we found 
more efficient to select randomly the SU{2) subgroups among the ^ ^~ -* subgroups acting 
on each 2x2 submatrix. At each site the SU{2) subgroup identified by the indices i,j 
(1 < i < j < A^) was updated with a probability P = -^^, so that the average number 
of SU{2) updatings per SU{N) site variable was n = aN. In our simulations we always 
chose a ^ 1, decreasing a when increasing A^. At A^ = 21, 30 we used a = 0.5. (We should 
say that our choices of the values of a came from a rough study of the performances of 
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the algorithm since the construction of the "most" efficient algorithm was not among our 
principal purposes.) Such an algorithm turns out to be quite efficient in the range of N and 
(3 (and correlation lengths) we considered. 

An important class of observables of the SU{N) x SU{N) chiral models can be con- 
structed from the correlation function G{x — y) = { ;^Re Tr [U{x)U{yy] ). 

The inverse mass gap C,w is extracted from the long distance behavior of the zero space 
momentum (wall- wall) correlation function constructed with G{x). Moreover we measured 
the diagonal wall-wall correlation length ^^ to test rotation invariance. M = 1/^^ should 
reproduce in the continuum limit the mass of the fundamental state. Another definition 
of correlation length ^g comes from the second moment of G{x). In the small momentum 
regime, p'^C.c <^ 1, we expect the behavior 

where G{p) is the Fouries transform of G{x), Mq = ^q^ and Zq is a constant. On the lattice 
we can use the two lowest components of G{p) to obtain the following definition of ^q'- 



e 



2_ 1 

G 



4 sin^ vr/L 



G(0,0) 
G(0, 1) 



(157) 



where G{kx, ky) is the lattice Fourier transform of G{x). 

In Table | we present data for the energy density E, the specific heat C = -^^, the 
magnetic susceptibility x — G{0), the correlation lengths ^g and ^^, the dimensionless 
ratios ic/^w and ^d/^w Data analysis were performed by using the Jackknife method. For 
comparison, in Table ^ we report some strong coupling results at A^ = oo obtained in 
Sections || and |Ti. 



Finite size effects were carefully checked (see also Ref. [Q). Finite size systematic errors 
in evaluating infinite volume quantities should be smaller than the statistical errors of all 
numerical results presented in this paper. In all cases we used lattice sizes L ^ 8.5^g- 

Figs. |l|, and ^ show the Monte Carlo data of the energy density and the specific heat 
respectively at A^ = 6, A^ = 9 and A^ = 15, with the corresponding strong and weak 
coupling series calculated in the previous sections. As in other asymptotically free models, 
at all values of A^ the specific heat shows a peak, connecting the two different asymptotic 
behaviors: monotonically increasing in the strong coupling region and decreasing at large /3. 
The position of the peak of G should give an estimate of the strong coupling convergence 
radius. 

As already observed in Ref. [Q and conffimed by our more recent data at A^ > 15, 
increasing A^ the peak of G moves slightly towards higher j3 values and becomes more and 
more pronounced, not showing any apparent convergence to a finite value. This might be 
an indication of a phase transition at A^ = oo, which a rough extrapolation would place at 
Psing ~ 0.305, with an uncertainty of few per mille. From A^ = 6 up to A^ = 30, we found 
the position of the peak to be very stable with respect to the correlation length: it occurs 
at Cg - 2.80 for A^ > 6. 

Rotation invariance at distances d ^ ^ is checked by measuring the ratio ^d/^w In Fig- § 
the ratio C,d/C,w is plotted versus C,g- Data show effective rotation invariance already for 
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^G — 2. At smaller correlation lengths, discrepancies from one are well reproduced by the 
Gaussian model (see Appendix B), which predicts 



6 iH* + \A+^ 



^'^ v^ln(-^ + ,/rT 



2V2^G ^ V 8?, 



•^a 



(158) 



Other important tests of scaling are based on the stability of dimensionless physical 
quantities, like the ratio ^g/Cw, which is in general different from one. Data of ^g/^w are 
shown in Fig. |^. Within statistical errors of few per mille, ^g/^w is stable for C,g ~ 2, and 
independent of N for A^ > 6, showing a rapid convergence to the N = oo value. We conclude 
that in the large A^ limit 






0.991 (159) 




with an uncertainty of about one per mille. As shown in Fig. ^, the comparison with the 
Gaussian model prediction 

(160) 

is very good at small correlation lengths .^g ~ 2, but then discrepancies arise since in the 
Gaussian case the continuum limit is one. 

Notice that scaling is observed even around the peak of the specific heat, even though its 
behavior with respect to A^ suggests the existence of a phase transition at A^ = oo. In a sense, 
at A^ = oo the modes responsible for the phase transition should be effectively decoupled 
from those determining the physical continuum limit, meant as the renormalization group 
trajectories where dimensionless quantities are stable. 

Similar considerations hold for the approach to asymptotic scaling, when using a "good" 
definition of temperature. A "good" definition of temperature turns out to he Te oc E (see 
Eq. ( |143| )), whose corresponding specific heat is, by definition, constant. From the Monte 



Carlo data of M = 1/^ui and the exact result (y), the effective A-parameters Al{N,(3) and 
A£;(A^, Pe) can be extracted. Fig. g shows the ratio Al{N, P)/Al,2i{N, /5), where Al,2i{N, (3) 
is the corresponding two loop function: AL^2iiN,/3) = {8tt/3) ' e"^""^. Around the region 
where the specific heat has a peak, also Al{N, P)/AL^2i{Af, P) shows a peak, whose shape 
tends to be singular when A^ — > oo for Psing — 0.305, consistently with the extrapolation 
from the specific heat. The observed peaks in the ratio Al{N, P)/Al^2i{N, P) correspond to 
dips in the /^-functions, and the singularity of Al(oo,/3) would reveal a singularity of the 
N = oo /3-function. The similarity with the behavior of the specific heat suggests that the 
peak of C and the dip of the /3-function have the same origin. They are presumably related 
to complex /5-singularities of the partition function close to the real axis , that for N ^ oo 
should finally pinch the real axis. 

As already observed in other contexts (see for example p3|j3^j35[] ) , the approach to 



asymptotic scaling gets an impressive improvement in the Pe scheme, where the dip of the 
/3-function disappears at all values of A^. Fig. |^ shows the ratio Ae{N , Pe) / Ae^2i{N , Pe) i 
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where Ae{N, jSe) is the effective A-parameter of the (5e scheme, and Ae,2i{N, [3e) is the cor- 
responding two loop function. Unhke Al{N,j3), Ae{N,(3e) appears to approach a smooth 
function Ae{oo, Pe), that is well approximated by the two loop formula Ae{oo,Pe) — 
{StcPe) e~^^^'^ even around the peak of the specific heat, which is placed at Pe — 0.220 
for A^ > 6 (the peak position being much more stable in N when described by the variable 
Pe instead of /5). 

In Section ^ we found that in the standard and the energy schemes the linear corrections 
to the two loop lattice scale in Eqs. ( |136| ) and (|149|) are small and of the same order of 



magnitude, although of opposite sign. Therefore perturbative arguments do not explain the 
failure of the standard and the success of the Pe scheme with respect to achieving asymptotic 
scaling. We believe the flattening of the specific heat to be the key feature of the Pe scheme. 
A coupling transformation eliminating the peak of the specific heat should move the complex 
/^-singularities away from the real axis, and therefore improve the approach to asymptotic 
scaling. 

The above considerations are consistent with the relationship (|146D . Indeed, assuming 
that at A^ = oo the /^-function of the Pe scheme is not singular (as data seem to show) 
and the specific heat has a divergence at Psing, then Eq. (|146|) predicts a singularity in the 
N = oo /5-function of the standard scheme at Psing- 

Asymptotic scaling within the Pe scheme can be also checked by using the strong coupling 
estimates of E, ^^ and ^g- In Fig- H we show the results using the strong coupling series 
obtained in Sections and |T| ( E up to O {p^'^), ^^ up to O (/3^) and ^g up to O (P^) ) 
at A^ = oo, and for Pe < 0.220, which should be the approximate position of the specific 
heat singularity. The comparison with the predicted mass-A parameter ratio is very good, 
especially for C,g whose series has alternate signs. 

Another quantity, we found to be very sensitive to the change of variable T -^ Te, 
is the renormalization constant Zg introduced in Eq. (|156|) . Its dependence on P can be 
determined by renormalization group considerations, indeed it must satisfy Eq. ( |116| ). One 
finds 



Zg = c{N) N"/'/^" z{T)-^ = c{N) p"^^'^'' 
= c{N) N^^/'"' ze{Te)-^ = c{N)Pe^'''" 



1 + O 




;i6i) 



where the function z{T) and ze(Te) were defined in Eqs. ( |140|) and ( |154|) , and c{N) is a 
constant independent of the regularization scheme. In Figs. P and |10| we plotted respectively 
the quantities cl{N,P) = ZgP"'^^''" and ce{N,Pe) = ZgPe ■ Fig. | shows also the strong 
coupling estimate of cl(oo, P) obtained from Eq. (|8^). cl{N, P) appears to be far from being 
approximately constant. It tends to be singular when A^ — > oo at Psing — 0.305, consistently 
with the previous estimates of the singularity. On the contrary, ce{N,Pe) shows a much 
smoother behavior and a rapid convergence at large A^. Again, perturbative arguments do 
not help to explain such different behaviors, since perturbative corrections are small for 
both schemes. Indeed from Eqs. (|140|) and (|154|) we can derive the following relationships 
at A^ = oo: 
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cl{oo,(3) = c(oo) 



1 + 



0.08521 0.01314 



P 



P' 



+ O 




(162) 



and 



ce{oo,Pe) = c(oo) 



1 + 



0.02271 0.00079 



Pl 



+ 



PI 




(163) 



We could get an estimate of c(cxd) from the data at the largest available Pe values, obtaining 
c(oo) ~ 0.10 (the perturbative corrections in Eq. ( |163|) are about 10% at (3e — 0.25). 

In Section |IIT| we found that at large A^ the correlation function breaks the Gaussian 
form at relatively high orders in the strong coupling expansion. Therefore in the very strong 
coupling domain, G{k) should be well approximated by Eq. ( 156 ) everywhere and not only 



for k"^ <C Mq. When the Gaussian approximation works we should find 



M{k) 



G{k) 
W) 



(mI + k' 



(164) 



Fig. |TT] shows the components (/c, 0) and (/c, k) oi M{ki, k2) at A^ = 30 and for various values 
of P chosen inside the expected convergence region of the strong coupling expansion. As 
expected from the strong couphng considerations, at /3 = 0.20 the propagator is effectively 
Gaussian. Discrepancies, although small, are clear for (3 > 0.28, especially in the diagonal 
components. 

In order to study the continuum limit of G{k), we consider the dimensionless function 



L(M,iv).a^^^-i. 

^ ' G{k-f3,N) 

L{k; P, N) is renormalization group invariant and therefore for k'^a'^ ^ 1 

L{k;P,N) ^ Ueec.N) , 



(165) 



(166) 



the function L(?/, A^) being independent of p. Notice that in the Gaussian model we 
would have L(?/) = y. In Figs. |1^ and |1^ we plotted the components (/c, 0) and (A;, k) 
of L{ki,k2; P,N) versus fc^^^ obtained respectively a.t N = 6 ior (3 = 0.31 and P = 0.33, 
and at A^ = 9 for j3 = 0.31 and j3 = 0.32. All sets of data follow a single curve for 
k'^^G ~ Qi/^'^1 which must be the continuum function L[y,N) at the given A^. Then dis- 
crepancies, i.e. scaling violations, arise. Comparing data at different P, we also learn that: 
(i) When there is rotation invariance at a scale p, i.e. L{p, 0; P, N) ~ L(p/\/2, p/\/2; (3, N), 
then L{p;P,N) ~ LijP'Q.^N). (ii) The diagonal components L{k,k; P,N) are closer to the 
continuum limit than the L{k, 0; P, N) ones. 

Fig. ^ gathers data of L{k, k; P, N) at different values of A^ and p. At relatively low 
values of k'^^Q, all sets of data follow approximatively a single curve, showing that L{y, N) 
rapidly converges to its large A^ limit. A^ = 6, 9 data give already a good approximation of 
L{y] oo). At larger k'^^Q, the scattering of the curves for different A^ is essentially due to 
scaling violations, because we used data at different ^g- 
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The extrapolation of L(y, A^) to negative y should provide another estimate of the ratio 
M/Mg = C,g/^w The solution yo of the mass gap equation L{y) + 1 = is |/o = —M'^/Mq. 
At sufficiently small y, L{y, oo) shows an approximate Gaussian behavior, indeed for y '^ 5 
it is well reproduced by a polynomial function y + by"^ with b very small: b ~ 0.01. By using 
this polynomial function to extrapolate, we find yo — — 1 within 1%, consistently with the 
more precise estimate (|159|) . 

We finally compare the Monte Carlo data of G{k) with the perturbative solutions G{k)p 
of the renormalization group equations (|138| ) and (|153| ). In Figs. [1^ and ^ we show 
data respectively ai N = 6 {P = 0.33) and N = 9 {P = 0.32), with their corresponding 
renormalization group equation solutions at the lowest order and at the next two orders, 
calculated in the standard and the Pe schemes (see Section ITVl ). In particular the lowest 
order solution of the renormalization group equation for G{k) is 

in the standard scheme, and 

G(*W=(^^(fln^] (168) 

in the Pe scheme. In G{k)p the scale is provided by the relationship 

where R is the mass-A parameter ratio in the corresponding scheme, and M = 1/^w is taken 
from the Monte Carlo simulation. 

Again, the best result comes from the Pe scheme, and in particular from its lowest 
order approximation. The A^; parameter shows a similar behavior too, its comparison with 
the corresponding two loop formula is already very good, and the O [Pe^] correction does 
not improve it. This seems to be a general feature of the Pe scheme: its lowest order 
renormalization group estimates do not need O (P'^^j corrections, since they are already 
very good. On the other hand, we should take into account that this scheme comes from a 
non perturbative change of variable. 

According to the conjectured S'-matrix ll|J^,|] and the large A^ factorization, 2-d SU{N) 
chiral theories should describe free particles in the limit N ^ oo. But when considering the 
two point Green's function G{x), we see that the realization of such physical properties is 
not trivial at all. Indeed, renormalization group considerations (see Section [TV] ) tell us that 
in the limit N ^ oo the asymptotic behavior of G{x) for x <^ 1/M is 

G(x) ^ W (-L) , (170) 

while a free Gaussian Green's function behaves as ln(l/x). Then G{x), at small x, seems 
to describe the propagation of a composite object formed by two elementary Gaussian exci- 
tations, suggesting an interesting hadronization picture: 

In the large N limit, the Lagrangian fields U, playing the role of non- interacting hadrons, 
are constituted by two confined particles, which appear free in the large momentum limit, 
due to asymptotic freedom. 
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APPENDIX A 

For low-dimensional representations the Z(^r) are (known) monomials in (3, while their 
counterparts Z(r+iN) may be reconstructed starting from the quantities 

Ti,,m, = {detuflTiU'^ nTrt/t'"^) , (Al) 

which in turn are recursively constructed by applying differential operators and algebraic 
Schwinger-Dyson equations (Ref. |TB[, Appendix D) to Fq 



Fo = ^(liv) = M2NP) + 0(/?3A^+2) . (A2) 

We have found the general form of a few interesting classes of character coefficients: 

d(ii-o)Z(ii.iN,^ ~ — - — Jn // _ i\| "^^+1 ' (^^) 

fc=0 [I' - i^)- 

{N(3y {N/jy-'^ 
dio;i)Z(o-i+iN) ^ —^ — Jn + -jizrw'^^^^ ■ '^^^^ 

Moreover we have explicitly computed the SU{N) character coefficients for the first few 
low-dimensionality representations 

rf(2;0)^(2;0) - liNpy + Jn+2 + 2N/3J'j, , (A7) 

d^i2,o)Zii2.o) ^ liNpy + Jn-2 + 2NI3J'^ , (A8) 

c^(i;i)^(i;i) ^ {Nfif - 2Jn + AN(3J'^ , (A9) 

dimHm - \iNl3f + \{NpfJM-i - \{NpfjM+i + NI3JN+2 - Jjv+3 , (AlO) 

rf(2,l;0)^(2,l;0) - \{N I3f + iV/3JjV-2 + [N^fj^^^ - {N^YJ^+i + NfiJ^^^ , (All) 
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rf(i3,o)^(i3;o) ^ l^Npf + J^_3 + NPJ^_, + l(iV/3)V;v-i - i(iV/?)V;v+i , (A12) 



C^(2;l)^(2;l) ^ ^iN f3f + ^(iV/?)V;V-l - 2iV/3J;v + 



1 - liNpr 



Jn+i — N(5Jn+2 



(A13) 



d(i2.i)Zi^2.i) ^ -{N(3f + N(3Jn-2 



1 - ^{N(3Y 



Jm-1 - 2NI3Jn - -{N(3fJr,+, 



(A14) 



APPENDIX B 



In order to improve our understanding, we briefly considered the Gaussian model, which 
bears a wide resemblance to the very strong coupling behavior of the chiral models for 
sufficiently large N. The Gaussian model is exactly solvable, and one finds that, for arbitrary 
values of the coupling and arbitrary direction 6 in the spatial lattice, Eq. (^) becomes 



and is solved by 



where 



p' + m' = , p^ = 4Y.^^^' (y) 



fi{6) = cos 6 SiYsh. [u cos 6) + sin6'arsh (z/sin( 



(Bl) 



[B2) 



m^ 



iy{e) = mW2 + — 



1 + 



\ 



1 — m'^{8 + m^) 



^cos2g ^ 

4 + m? 



-1/2 



(B3) 



Introducing the auxiliary strong coupling variable 

1 

7 



4 + m2 
(related to /3 by /3 = 7 + 0(7^), odd function of 7), we obtain 

1 VI - 167^ 



m 
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1 + Jsm\2e) + 1672cos2(20) 



(B4) 



(B5) 



that is, z/ is an odd function of 7 for all 6 y^ 0. The crucial observation however concerns 
the limit 6^0, where 



KO) = ^\/l-47 



(B6) 
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which is not a function of 7 with definite parity. Therefore, while for all 6' 7^ the quantity 

p{e) = fi{e) + (cos^ + sin^)ln/3(7) (B7) 

is an even function of 7, this property does not hold at 6' = 0. When we consider Green's 
functions, we recognize that in strong coupling 

G{xi,X2) = /3"^+"^/7(xi,X2;/5) (B8) 

where H is always an even function of /? (and 7). As a consequence, the purely kineniatical 
singularity at ^ = in Eq. (p5| ) prevents the general relationship 

p(e) = - lim ^^ B9 

from holding at 6* = 0. For all other orientations exponentiation can be shown to hold, and 
in particular the principal diagonal correlation (a quantity of fundamental relevance for tests 
of rotation invariance) has the proper exponential decay and the corresponding mass gap pd 
can be shown to coincide with the value extracted from the diagonal wall-wall correlation: 

^, = - lim '■'G^...fe-^^ = l^l/V2) . (BIO) 

\x\^oo \x\ 

Notice that for arbitrary 6 the two definitions do not strictly coincide. 



APPENDIX C 

We could establish the following recursive relationships 

Co{xi,X2) - Co(2;i - 1,0:2) - Co{xi,X2-l) = , (CI) 

with boundary conditions 

Co(xi,0) = Co(0,X2) = 1 (C2) 

and (well known) solution 

Co(x„X2) = l^^'l^') - [^'t^') ; ^^3) 

C2(X1,X2) - C2{Xi - 1,^2) - C2{Xi,X2 - 1) = 

Co{Xi - 1, X2 + 1) + Co{Xi + 1, X2 - 1) - Co{Xi - 1, X2) - Co{Xi,X2 - 1) , (C4) 

with boundary conditions 
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and solution 



C2{Xi,X2] 



C2{XI,0) = 


a;i(a;i + 1) , 


^2(0,^2) = 


3:2(^2 + 1) , 


' xi + a;2 \ 


a;i(a;i + l) ^ 0:2 (3^2 + 1) 
_ 0:2 + 1 Xi + 1 _ 


. ^1 / 



(C5) 
(C6) 



C4{Xi,X2) - C4{Xi-l,X2) - C4{Xi,X2-l) = 

C2{xi + 1, a;2 - 1) + 6*2(3:1 - 1, 0:2 + 1) - 6*2(3:1, X2) + Co{xi + 2, 0:2 - 1^ 
+ Co{xi - 1, a;2 + 2) + Co{xi, 2:2 - 2) + Co{xi - 2, X2) - 3Co(xi, 0:2) , 

with boundary conditions 

^.(.1,0) = "'^^'^' +^^^^%^ + '"^' 



(C7) 



(CJ 



and solution 



C4(Xi,X2) 



Xi + X2 
Xi 



X1X2 + 2xi + 2x2 



a;2(a;2 + 3) a;i(a;i + 3) 



xi + 1 a;2 + 1 

(a;2 - l)x2(a;2 + l)(a;2 + 2) (xi - l)a;i(a;i + l)(a;i + 2) 2xiX2 



+ 



2(xi + l)(xi + 2) 2(x2 + l)(x2 + 2) 

The function B can be computed by the relationship 

B{xi,X2) = ^aC,^"-'(xi,X2) , 



Xi + X2 



■ (C9) 



(CIO) 



.(«) 



where Cg (xi,X2) is the number of minimal self-avoiding paths connecting the origin with 
X and forming a right angles. Cq satisfies the normalization condition 

^Cq (xi,X2) = Co(Xi,X2) . 
a 

Solving appropriate recursive equations it is possible to prove that 

y(a)( \,xi ,xo f tl . h \ I tlt2 



(Cll) 



J2 C!f\x^,X2)trt 



X2 
2 



XI, X2 



and as a consequence 



1-ti ' i-t2y \{i-ti){i-t2 

g + l 

tit 



for even a , 



h 



for odd a 



J2 B{xi,X2)tlH. 



X2 
2 



2tit2 



11, 12 



'l-h-t2f 



We then trivially obtain 



-B(xi,X2) = ■ Co(Xi,X2) . 

Xi + X2 



(C12) 



(C13) 



(C14) 
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APPENDIX D 



We report here the N = oo strong couphng series of some relevant quantities: 
E = 1 - /5 - 2/5^ - 6/3^ - 38/3^ - 240/3^ - 1812/3^^ - 14126/3^^ + O (/3^^) , 



:di) 



X = 1 + 4/5 + 12/3^ + 36/5^ + 100/5^ + 284/5^ + 796/3^ + 2276/5^ 
+6444^^ + 18572/3^ + 53284/?^° + O (f3^^' 



:D2) 



Zg = P'^\l-P^- 6/3^ - 4/3^ - 20/3^ - 24/3^ - 148/5^ - 240/5^ + 0(p 



?io 



(D3) 



M^ = /3-1 [l - 4/3 + 3/3^ - 2/3^ - 4/3^ - 12/3^ - 40/3^ + 84/3^ - 320/5^ + O (/?^°)] , (D4) 



^is^de = - In /5 - 2/5 - ^/^^ - 2/5^ - ^/^^ - 8/5^ - ^/^^ - 84/^8 ^ ^ (^^9^ p5) 



l^diag 



V2 



6 



(D6) 
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FIGURES 

FIG. 1. Energy and specific heat versus (3 for A^ = 6. The full and dashed lines represent 
respectively the strong coupling (up to O (/3^^) for E and up to O {(3^'^) for C) and the weak 
coupling series. 

FIG. 2. Energy and specific heat versus (3 for A^ = 9. The full and dashed lines represent 
respectively the strong coupling (up to O (/9^^) for E and up to O (/3^^) for C) and the weak 
coupling series. 

FIG. 3. Energy and specific heat versus /3 for N = 15. The full and dashed lines represent 
respectively the strong coupling (up to O {/3^'^) for E and up to O (/3^^) for C) and the weak 
coupling series. 



FIG. 4. The ratio ^d/^w versus ^g- The full line represents the Gaussian prediction (|15j 



FIG. 5. The ratio ^g/^w versus ^g- The full line represents the Gaussian prediction ( |16C| ). The 
dashed line is the result of a fit. 

FIG. 6. A.l{N , j3) / A.L^2i[N , (3) versus (3. The dotted lines connecting different sets of data are 
drawn to guide the eyes. 

FIG. 7. Ke{N,(3e)/^e,2i{N,Pe) versus /3e. 

FIG. 8. Asymptotic scaling test by using strong coupling estimates. The dotted line represents 
the exact result (^. 

FIG. 9. cl = Zg(3'^^' versus /?. The dashed line shows the strong coupling series for N = oo. 
The dotted lines connecting different sets of data are drawn to guide the eyes. 

FIG. 10. CE = ZgPe^^° versus Pe- 

FIG. 11. M{k, 0) and M{k, k) versus P at iV = 30. 

FIG. 12. L(A;,0) and L{k, k) versus fc^^^ at iV = 6 and for (3 = 0.31 and (3 = 0.33. The dotted 
lines connecting different sets of data are drawn to guide the eyes. The dashed line represents the 
Gaussian prediction. 
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FIG. 13. L{k,0) and L{k, k) versus k'^5,G at TV = 9 and for (3 = 0.31 and (3 = 0.32. The dotted 
lines connecting different sets of data are drawn to guide the eyes. The dashed Hue represents the 
Gaussian prediction. 



FIG. 14. L{k,k) versus k^if.. 

FIG. 15. G{k, 0) and G{k, k) versus /c^^g at iV = 6 and /? = 0.33. The low and high sets of lines 
(full, dashed and dotted) show respectively the renormalization group predictions of the standard 
and the (3e schemes. Full, dashed and dotted lines represent respectively the lowest, and the one 
and two next order perturbative solutions. 

FIG. 16. G{k, 0) and G{k, k) versus k'^^l. at iV = 9 and /? = 0.32. The low and high sets of lines 
(full, dashed and dotted) show respectively the renormalization group predictions of the standard 
and the (3e schemes. Full, dashed and dotted lines represent respectively the lowest, and the one 
and two next order perturbative solutions. 
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TABLES 



TABLE L Summary of the numerical results. Errors of data marked by an asterisk could be 
underestimated. When more than one lattice size appear, the corresponding results were obtained 
collecting data of simulations at the reported lattice sizes (which were, in all cases, in agreement 



within the errors). 
















/3 N L 


E 


C 


X 


?G 


(.n, 


ia/in, 


Cd/5» 


0.20 15 18,21 
0.20 21 18,21 
0.20 30 18 


0.781405(7) 
0.781427(7) 
0.781422(10) 


0.0527(3) 
0.0527(3) 
0.0522(6) 


2.9380(7) 
2.9381(5) 
2.9379(8) 


0.786(3) 
0.788(2) 
0.787(2) 


0.8360(4) 
0.8358(3) 
0.834(2) 


0.941(3) 
0.943(2) 
0.944(3) 


0.9724(4) 
0.9726(4) 
0.973(2) 



0.28 
0.28 
0.28 



15 
21 
30 



18,24 

18,24 

18 



0.65000(5) 
0.65290(3) 
0.65352(3) 



0.191(5) 
0.170(3) 
0.163(4) 



7.249(6) 
7.069(5) 
7.032(5) 



1.560(4) 
1.532(3) 
1.529(3) 



1.587(2) 

1.5605(10) 

1.5544(15) 



0.9834(13) 
0.9819(11) 
0.9837(12) 



0.9922(6) 
0.9917(4) 
0.9916(5) 



0.29 
0.29 



9 
30 



24,30 
21 



0.58772(7) 
0.63058(3) 



0.412(5) 
0.208(6) 



13.32(3) 
8.497(8) 



2.369(10) 
1.740(4) 



2.395(9) 

1.765(5) 



0.9889(9) 
0.9857(7) 



1.000(2 
0.992(2 



0.295 



9 24,30,36,42 0.56278(6) 0.444(5) 



18.03(4) 



2.913(9) 2.949(13) 0.988(2) 



0.996(3 



0.30 
0.30 
0.30 



15 
21 
30 



24, 



30,36,42 

24 
24 



0.56806(4) 

0.58799(14) 

0.59927(8) 



0.68(2) 
0.65(3) 
0.38(2) 



16.574(13) 

12.91(3) 

11.35(2) 



2.742(5) 
2.309(7) 
2.114(7) 



2.767(9) 
2.333(7) 
2.137(10) 



0.9907(10) 

0.9899(7) 

0.989(2) 



1.000(1 
0.997(2 
0.993(3 



0.3025 
0.3025 



21 
30 



30 
30 



0.56525(19) 
0.58479(10) 



1.02(5)* 
0.79(5)* 



17.02(5) 
13.24(3) 



2.786(10) 
2.338(7) 



2.813(12) 
2.362(9) 



0.9903(8) 
0.990(2) 



0.997(2 
0.999(2 



0.304 



30 



24 



0.5827(4) 



2.6(3)* 



17.52(9) 



2.839(8) 2.866(11) 0.9907(13) 0.998(2 



0.305 
0.305 
0.305 



15 
21 
30 



36 
30 
30 



0.53415(8) 

0.54098(10) 

0.54658(14) 



0.52(2) 
0.73(3) 
1.05(10)* 



26.86(7) 
24.14(6) 
22.13(5) 



3.782(16) 
3.523(10) 
3.320(12) 



3.82(3) 
3.57(2) 
3.35(2) 



0.989(3) 
0.988(2) 
0.990(2) 



0.996(4 
0.998(3 
1.000(3 



0.31 
0.31 
0.31 
0.31 



6 

9 

15 

21 



60 

48 
45 
42 



0.48187(2) 
0.50030(4) 
0.51178(4) 
0.51548(6) 



0.257(3) 
0.302(6) 
0.354(7) 
0.41(2) 



65.2(2) 
47.25(12) 
39.06(10) 
36.66(12) 



6.63(3) 

5.44(3) 
4.80(2) 
4.61(2) 



6.69(4) 
5.49(5) 
4.84(3) 
4.65(3) 



0.9909(10) 
0.9908(17) 
0.9911(12) 
0.9915(10) 



1.003(3 
0.998(5 
1.001(2 
0.999(3 



0.32 
0.32 



9 
15 



66 
66 



0.47234(3) 
0.48072(2) 



0.252(3) 
0.264(5) 



82.4(3) 
70.5(3) 



7.70(4) 
6.96(4) 



7.78(7) 
7.04(4) 



0.990(2) 
0.9892(11) 



1.002(3 
1.000(3 



0.33 



102 



0.43706(2) 0.212(3) 175.3(2.0) 12.13(16) 12.23(21) 



0.991(3) 1.001(9 



TABLE IL Strong coupling results at N 



oo. 



/3 


E 


C 


X 


ia 


£,n, 


0.20 


0.7814220 


0.052534 


2.93042 


0.78783 


0.83547 


0.25 


0.7090132 


0.101253 


4.53371 


1.14122 


1.17106 


0.28 


0.6556834 


0.154155 


6.24293 


1.52257 


1.51843 


0.29 


0.6352185 


0.179292 


7.01904 


1.71989 


1.68452 


0.30 


0.6129234 


0.210053 


7.93314 


1.99083 


1.89275 
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